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1. INTRODUCTION 

The current modern electrical power system is one of the most complex systems in the human 
history. Classical generation power stations are centralized, monitored, built away from human cities and 
transmitted over thousands of kilometers to the required electrical loads. In the recent decades the electrical 
power system suffered from different blackouts due to different reasons which affected millions of customers 
all over the world. To overcome some of modern grid problems which includes but not limited to power 
losses, outages, blackouts, voltage flickers, harmonics, efficiency, voltage instability and costs decentralized 
distributed generation (DG) sources are used close to the customers connected or isolated from main 
electrical grid which provides another solution of providing less polluted electrical energy sources as most of 
the DG technologies are renewable environmental friendly sources such as photovoltaic systems and wind 
turbines [1-7]. Another reason for involving DG technologies in some countries is to make a diversity in the 
companies supplying electricity, thus establishing a competitive electrical market where the customer 
benefits from different categories and prices of electrical bill. The Department of Energy (DOE) considers 
the DG as the generation of electricity from the sources which are varying from less than a kilowatt (kW) to 
tens of megawatts (MW) in size for single distributed generation unit [8]. The merits of placing and installing 
of distributed generation (DG) units in radial distribution network beside the economic and environmental 
benefits includes voltage profile enhancement, system total power loss reduction, power quality 
improvement, strengthen system reliability, protection, increase system stability, system better performance 
and power flow control [8-9]. It is worth mentioning that the improper sizing and random locating of DG 
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units in the power system can lead to increase the active power losses in the distribution system. 
DG resources can be categorized as mentioned in [10] into four types; 

DG type I: capable of injecting active power in kW and operates at unity power factor. Examples of this type 
are photovoltaic, fuel cells and micro-turbines. 

DG type II: capable of injecting reactive power in kVAR. Reactive power compensators, capacitors and 
synchronous compensators belong to this type. 

DG type III: capable of injecting both active and reactive power in kW and kVAR respectively. DG 
resources depending on synchronous machine as cogeneration are part of this type of distributed generation. 
DG type IV: capable of injecting active power in kW and consuming reactive power in kVAR at the same 
time. The obvious example of this type of DG is induction generators used in wind farms. 

Many attempts have been done using artificial intelligence (AI) based techniques to optimally size 
and locate the Distributed Generation (DG) in the distribution systems. The whale optimization algorithm 
and its modifications belong to the swarm intelligence methods which are a part of the AI based methods. 
Most swarm intelligence techniques are inspired by the creatures and species in nature. 

The application of ant lion optimization algorithm to optimally allocate and size renewable DG 
technologies (wind turbines and photovoltaic systems) are discussed in [11] where the loss sensitivity factors 
are used to determine the required candidate buses to connect the DG units. However, it can be noticed that 
when installing two wind turbines on 33 bus test system to minimize the power loss in kW, 
the sizes of the DG units used in KVA are bigger than the ones used in [12] on the same test system 
by 27.7%. 

Ref. [13] proposes a meta-heuristic algorithm called invasive weed optimization algorithm to find 
the optimal sizing of the placed DG and location of the bus where the DG would be installed is determined 
by loss sensitivity factor. The issue in this study that the reduced total real power loss after placement of the 
DG is not the optimum compared to other methods applied to the same test systems. 

Ref. [14] focused on placement of different DG renewable resources such as wind turbines and PV 
arrays in radial distribution feeders, the study used type I, type II and type III DGs. The optimization method 
used to reach the optimal global solution is the whale optimization algorithm WOA with the assistance of 
power loss index. The introduced technique is tested on small and large test feeders 15 bus, 33 bus, 69 bus, 
85 bus and 118 bus test systems. The disadvantage of this work is that the results obtained from testing the 15 
bus and 118 bus test systems are not compared with other techniques and methods previously applied to the 
same test feeders, so the effectiveness of the technique used is not verified. Regarding the 33 bus and 69 bus 
test systems the technique couldn’t reach the optimal voltage profile improvement as for the instance an older 
method proposed in [15] reached better system voltage profile when applied to the same test feeders. 

Ref. [16] introduced a modern optimization technique to approach the solution of the DG allocation 
and sizing problem in distribution feeders. The stud krill herd algorithm is utilized to get the minimum power 
losses of the system under study. The algorithm is implemented and tested on 33 bus, 69 bus test systems and 
the 94 bus Portuguese radial distribution system. The study placed one DG, two DG and three DG to prove 
the abilities of the technique used. The results show that the algorithm achieved its target to reach lowest 
possible active power losses in the system at the expense of using larger sizes of placed DG units when 
compared to other optimization algorithms in references [17-19] applied to the same 33 bus, 69 bus test 
systems and the 94 bus Portuguese radial distribution feeder studied in the last previous years. 

Ref. [20] presented the grey wolf optimizer as a method to solve the problem of optimal locating 
and sizing of multiple DG units in radial distribution system. The method is implemented on IEEE 69 bus 
test system to verify its effectiveness and capabilities. The system power losses were enormously reduced 
with the placement of different DG sizes accompanied by enhancement in the system voltage profile. 
However, the average elapsed time to converge the optimal solution is higher than other methods applied to 
the same test system as the bat algorithm and gravitational search algorithm. 

In this paper, the modified whale optimization algorithm is applied to three different IEEE radial 
distribution feeders. This proposed technique is intended to find the optimal solutions to allocate and size DG 
units (type I, type II and type IID) in radial distribution feeders. The algorithm follows a set of predetermined 
steps to find the global optimal solution to this kind of hard optimization problems. This proposed method is 
implemented for the first time to solve this kind of optimization problems which includes sizing and locating 
DG units in any radial distribution feeder. This novel algorithm has the capability to find the optimal solution 
regarding the lowest system minimum real power losses accompanied by voltage profile improvement 
without violating the problem constraints such as power conservation boundaries, bus voltage constraints, 
DG real and reactive power limits and maximum thermal capacity of branch lines constraint. 

The rest of the paper is divided as follows: Section 2 presents the traditional whale optimization 
algorithm and the modifications applied to it to be modified and enhanced. Also, the solution algorithm to 
implement the modified whale optimization algorithm to the problem of optimal allocation and sizing of 
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different DG units in radial distribution feeder is explained in detail. Section 3 presents the implementation 
and computer based program simulation results compared to other optimization techniques using MATLAB 
facilities. Finally Section 4 includes the conclusions. 


2. RESEARCH METHOD 
2.1. Whale Optimization Algorithm (WOA) 

In nature, the real whales spread all over the world in oceans and seas. Scientists studied their way 
of living and hunting. In Ref. [21, 22], the traditional WOA is a nature inspired meta-heuristic optimization 
algorithm which mimics the social behaviour of humpback whales. The bubble net hunting method used by 
humpback whales in hunting its prey such as group of fishes close the surface of the water is utilized by the 
artificial whales in the proposed algorithm. The method is simple and easy to be implemented. 
The mathematical model could be described as follow: 

A. Search for prey 

The first step is exploration phase, in this phase the humpback whales search randomly for the prey, 
so the variation of vector A is used to represent this searching process. The position update of each whale is 
done by moving towards a random known search agent, so the global search happens. As shown is (1) and 
(2) describe the previous behavior: 


D=IC. Xrand - X| (1) 

X(t+1) = Xrand - A.D (2) 
where 

A= 2.ar—a (3) 

Calr (4) 


A is coefficient vectors 

C is coefficient vectors 

t is the current iteration 

X is the position vector 

a is linear decreasing from 2 to 0 

r is the random number between [0,1] 


B. Encircling the prey 

The location of the prey is recognized by humpback whales which encircle them. The other search 
agents in the search space change and update their positions towards the defined best search agent. 
The behavior is described mathematically as follows: 


D=| C.X*(t) —X(t)| (5) 
X(t+1)=X*(t) — A.D (6) 


where : 
X* is the best value of the position vector 


C. Bubble-net attacking method 

To hunt the prey the humpback whales makes two mechanisms simultaneously which are shrinking 
circle and moving in spiral path shape to update their position. After determining the distance between the 
whale and prey, the movement takes place as shown in Figure 1. An assumption is made that the probability 
to choose between either to shrink encircle or to move in spiral shaped path is 50%, so a variable (p) is 
introduced to indicate that probability, Another constant is introduced to represent the spiral shape path called 
(b). The behavior is described by the following equation: 


X*@)=A.D iffp<05 


sea {ped cos(2.7.1) + X*(t) iff p>0.5 


(7) 
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where 
(p) is arandom number between [0,1]. 
(1) is random number from [-1,1]. 





Figure |. Bubble-net feeding behaviour of humpback whales 


2.2. Modified Whale Optimization Algorithm (MWOA) 

In Ref. [23], procedures are presented in order to modify and enhance the traditional whale 
optimization algorithm regarding the performance and results obtained. An inertia weight (w) € [0, 1] is 
added to introduce the modified whale optimization algorithm (MWOA), in this case (5) and (6) are modified 
as follows: 


D =|C.@X* (t )- X(t)| (8) 
X(t +1) =0X* (t )- A.D (9) 


Where: 
(@ ) an inertia weight €[0,1] 


Also (7) representing the update of the position in the spiral path would be modified as follow: 


wX* (t) —A.D iff p<0.5 


Xt+)= i oes cos(2. 7.1) + wX* (t) iff p>0.5 


(10) 


The modified algorithm steps of the MWOA can be summarized as: 

a) Initialize the whales’ population Xi where i = (1, 2, 3, ., n) and set the Maxgen (maximum number of 
iterations). Set t = 1. 

b) Calculate the fitness value of Xi where i = (1, 2, 3, ...., n), identify the best search candidate solution 
X*, 

c) Repeat the following: For every Xi (i = 1, 2, 3, ... , n), update p, C, A, 1, a. If p < 0.5, then if | A I< 1, 
update the position of the current search agent by (8) and if | A |= 1, select a random search solution 
Xrand and update the position of the current search agent by (2). But when p = 0.5, update the position 
of the current search by the (10) Check if any search agent goes beyond the known search and modify 
it. Calculate the fitness of Xi (i= 1, 2 , 3, ... , n) , and if there is a better solution, find the best search 
solution X*. Set t= t+1. Until t reaches maximum number of iterations, the algorithm is finished. 

d) Return the best optimization solution X* and the best optimization value of fitness values. 


2.3. The Power Flow Calculations 

From [7] and others, Newton Rapshon, fast decoupled load flow, backward forward sweep and 
Gauss Seidel are some of the most popular methods to solve the load flow problem of radial distribution 
feeder. In this study, the Newton Raphson technique is used to find the power flow solution for each test 
feeder. Figure 2 shows the single line diagram of a general radial distribution feeder model. The governing 
load flow equations used could be stated as follow [24]: 


Pint = Py — rin (Pi27+Qi7) / Vi? = Pritt (11) 
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Qin = QHxin (Pi? + Qi?) / Vi? = Quist + Qciv1(12) (12) 
Vier? =Vi2-2 (ris Pit Xin 1 Qi) +((tis 12+ Xin 17)*(Pi24+Q; 7) ) /Vi2 (13) 
Piossi1=Tin (Pi +Q?7)/ViP (14) 
Quossit=Xie1 (Pi? +Q?)/Vi? (15) 


Where, 

P;: active power flows into the sending of branch i+1 connecting node i and node i+1 
Qi: reactive power flows into the sending of branch i+1 connecting node i and node i+1 
Vi: bus voltage magnitude at node i 

Qci+i : reactive power injection from capacitor at node i+1 

PLis1: active power load connected to node i+1 

Qti+i: reactive power load connected to node i+1 

ri+1 ; line resistance between node i and node i+1 

Xi+1! line reactance between node i and node i+1 

Prossit1; active power losses of branch connecting node i and node i+1 

Quossi+1; feactive power losses of branch connecting node i and node i+1 


Vo k Pk, ok Wei Wn 


Po. Qo A ile Ft ri 
eet ee \e 








i 


Figure 2. Single line diagram of a model of radial distribution feeder 


2.4. Objective Function 

The purpose of optimal sitting and sizing of DG units in the radial distribution system is to enhance 
the voltage profile of the system and decrease the real power loss in kW on condition not to violate the 
predetermined system constraints. From [17] the optimal solution is to minimize the real power losses which 
can be stated as: 

Minimize F = min (Pyotal Joss) (16) 


where 
P.otal loss 1S the total real power losses of distribution network 
According to the following constraints: 
A. Voltage constraints 
Vmin<Vi< Vmax (17) 
where 
Vmin: minimum bus voltage limit 


Vmax: maximum bus voltage limit 


B. Power balance constraint 
Piasech Pe Passe Dyes (18) 
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Where, 

P3jack 1S the slack bus real power 
Poss is the system power loss in kw 
P, is the system load in kw 

Ppyg is the power of DG unit 


C. Upper and lower limits of DG unit 
Poa” SS Pog S Prem (19) 


P5q” minimum DG unit size 
DG Maximum DG unit size 


From [25], the upper and lower limits can be set as follows: 

60< Ppe< 3000 (20) 
Where the limits are in kVA, kVAr and KW for type III, type II and type I DG respectively 
D. Current constraints 

Ti < Jimax for i=1, 2,.... Nb (21) 


Ti represents the branch current of the i™ branch of the network consisting of Np» number of 
branches. Jimax is the maximum permissible current flows through the i branch. 


2.5. Algorithm Steps 

a) Using Newton Rapshon method solve the load flow problem for the given test feeder and determine the 
system total power loss. 

b) Initialize the random number of agents, set the iteration counter = 1, set the values for minimum and 
maximum size of each DG unit, specify the number of DG units used. 

c) Determine the real power los for generated population by performing load flow calculations. 

d) Select the DG value with the lowest power losses as the current best solution. 

e) Ifp<0.5 and A>= 1, Update the position of the agents using (1) and (2). 

f) Repeat step c. 

g) Replace the current best solution with the updated values if the obtained results are lower than the 
current best solution. Otherwise go to step h. 

h) If p< 0.5 and A < 1, Update the position of the agents using (8) and (9). 

i) Repeat step c. 

j) Replace the current best solution with the updated values if the obtained results are lower than the 
current best solution. Otherwise go to step k 

k) If p> =0.5, update the position of the agents using (10). 

1) Repeat step c. 

m) Replace the current best solution with the updated values if the obtained results are lower than the 
current best solution. Otherwise go to step n. 

n) Start the second iteration. 

0) Repeat until the maximum number of iterations is performed and print the results. 


3. RESULTS AND ANALYSIS 
3.1. The First Feeder 

The first test case is the 15 bus test feeder with the single line diagram shown in Figure 3, the radial 
distribution feeder data is given in [25]. The load flow calculations applied to that test feeder before adding 
any DG units to the system. The power factor is 0.7. The total apparent power is 1752 kVA. 
The reactive power losses and the total active power losses are 57.297 kVar and 61.954 kW respectively, 


the minimum and maximum buses voltages range is 0.946 SViF1.0 p.u. 

The modified whale optimization algorithm (MWOA) technique is implemented and the DG units 
have been sized and placed in the system. The results seem to be superior regarding the real power loss which 
has been decreased enormously and the system voltage profile has been improved. 
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Table 1 presents a comparison between the results of placing DG type I of different techniques. 
The superiority of the proposed method in this paper is clear. The technique results compared to ant lion 
optimization and traditional whale optimization algorithm showed that it produced the lowest real power 
losses with loss reduction percentage up to 38.68% which is presented in Figure 4 and the voltage profile is 
improved as the minimum voltage bus in pu in the system at bus number 13 is 0.966. 

Table 2 presents another comparison of results between placing different DG types. It is obvious 
that type III DG has the lowest system real power losses with reduction in power losses up to 64.85% and 
best enhanced voltage profile with minimum voltage bus 0.97 pu. 

From Table | and Figure 4 it can be deduced that the proposed method has a superior performance regarding 
decreasing the real power loss and improving system voltage profile compared to other optimization 
techniques. 






Substation 


Figure 3. Single line diagram of 15 bus test system 


3.2. The Second Feeder 

The second tested feeder using the proposed modified whale optimization algorithm MWOA is the 
33 bus system. It consists of main feeders and three laterals. The system single line diagram is shown in 
Figure 5. This test feeder has a total load of 3720 kW and 2300 kVAr at a voltage level of 12.66 kV. 
The system data are presented in Ref. [14]. The configuration of the system before installing any DG unit is 
as follow: the real power loss in kW is 210.9 and the minimum bus voltage is 0.9038. The results of using 
three different DG types (type I, type II and type III) applied to the test feeder are shown in Table 3. From 
Table 3, it is clear that DG type II installation results in the most reduction in real system power loss 
compared to the other DG types. Also the best minimum bus voltage for the tested feeder was achieved when 
installing DG type IIL. 

Table 4 presents a comparison between different methods and techniques applied to this test feeder, 
it can be deduced that the proposed method has the lowest active power loss in Kw and the most improved 
system voltage profile. Figure 6 shows the enhancement in system voltage profile after adding DG type III. 
Also in Figure 7 the real power loss has the lowest value when the proposed technique is applied to the test 
feeder and compared to other methods including the classical whale optimization algorithm. 
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Figure 4. Power loss reduction percentage of 15 bus test system 
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Table 1. Optimal Results of DG Type 1 Placement in 15-Bus Feeder 
After Single DG Placement 








B DeforeG aren 
Parameter Placement Optimization Whale Optimization Algorithm Proposed 
(ALO) (WOA) Method 
[27] [15] 
Year of technique publication June 2017 May 2017 May 2018 
(Reference) 
Final real power 
loss(kW) 61.95 45.803 42.818 37.988 
Final reactive power 
loss (KVAR) 57.29 41.88 37.97 34 
Loss Reduction % for real power losses 26% 30.88% 38.68% 
Minimum Bus Voltage 0.94 0.952 0.959 0.966 
DG size (KW) 675.24 673.86 948.29 
DG location 6 15 3 








Figure 5. Single line diagram of 33-bus test system 


Table 2. Optimal Results 15 Bus Feeder with DG Type 1, Type 2, and Type 3 
After Single DG Placement 








Parameter ¥ Dera DG Type I DG Type II Be Tye a 
Placement ee Sins Size in KVA 
Size in KW Size in KVAR 
Power factor 0.9 lag 
Final real power 
loss(kW) 61.95 37.988 37.221 21.773 
Final reactive power 
loss (KVAR) 57.29 34 33 18.29 
Loss Reduction % for real power losses 38.68% 39.92% 64.85% 
Minimum Bus Voltage 0.94 0.966 0.965 0.97 
DG size 948.29 963.4 1040.14 
DG location 3 3 4 





Table 3. Optimal Results 33 Bus Feeder with DG Type 1, Type 2, and Type 3 
After Single DG Placement 








Parameter Beene DG Type I DG Type II os Type a 
Placement ae ae Size in KVA 
Size in KW Size in KVAR 
Power factor 0.9 lag 
Final real power 
loss(kW) 210.9 103.94 143.588 64.34 
Final reactive power 
loss (KVAR) 143.03 75 96 50 
Loss Reduction % for real power losses 50.72 31.91 69.49 
Minimum Bus Voltage 0.9038 0.954 0.926 0.968 
DG size 2575.35 1252.53 3000 
DG location 6 30 6 
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Table 4. Optimal Results of DG Type 1, Placement in 33 Bus Feeder 
After Single DG Placement 
Before DG Backtracking Stud Krill Herd Whale Optimization 
Parameter tae Search Algorithm Algorithm Algorithm Proposed 
(BSOA) (SKHA) (WOA) Method 
[13] [17] [15] 
Year of technique publication 015 2017 May 2017 May 2018 
(Reference) 
Final real power 
loss(kW) 210.9 118.12 111.01 125.161 103.94 
Final reactive 
power loss 143.03 = 81.71 89.28 75 
(KVAR) 
Loss Reduction % for real power losses 44.01 47.36 40.65 50.72 
pees 0.9038 0.9441 0.9424 0.9272 0.954 
Voltage 
DG size (KW) 1857.5 2590 1542.67 2575 
DG location 8 6 30 6 
Voltage Profile Before and After Adding DG 
1.02 
1 
x ) 
2 0.98 “ j x 
7 2 hy Itage i 
= = voltage in pu 
en ¥ ' \ before adding DG 
5 = ! = 
> 0.94 \ r N —— voltage in pu 
nos \ ~——-+4 ~~ after adding DG 
~ 
0.9 
0 5 10 15 20 25 30 35 
Bus No. 
Figure 6. The effect of installing DG type III on system voltage profile 
Power loss 
in KW 
140 + 
120 
100 
80 : 2 E : : ——s 
60 
40 J : =k pL E u = = 
20 @ Power loss 
Oo T T : in KW 
Backtracking Stud krill Herd Whale Proposed 
Search Algorithm Optimization MWOA 
Algorithm (SKHA) Algorithm method 
(BSOA) [17] (WOA) 
[13] (15] 
Figure 7. Results of system power loss reduction in kW of 33- bus test system 
4. CONCLUSION 


In this paper, a modified whale optimization algorithm has been successfully implemented for 
optimal location and sizing of DG units type I, type II and type II in various distribution test systems. 
The designed problem has been formulated as an optimization task with computing system power losses and 
enhancing system voltage profile. The effectiveness of the proposed technique is presented by using different 
test systems. The results have been compared regarding using different DG types and also with those 
obtained using other optimization methods. It is clear from the comparison that the proposed optimization 
method provides a superior performance in terms of system real power losses and improving system 
voltage profile. 
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